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I. INTRODUCTION 



Loop Quantum Gravity (LQG) has become in the past years a mayor player as a candidate for 
a quantum theory of gravity [3, 0] • On the one hand it has matured into a serious contender 
together with other approaches such as String/M Theory, but on the other is it not as well 
understood, neither properly credited as a real physical theory of quantum gravity The 
purpose of this contribution is twofold. First, we would like to provide a starting point for 
those interested in learning the basics of the theory and to provide at the same time an 
introduction to one of its application, namely the treatment of black holes and its entropy 
counting. 

Loop quantum gravity is based and two basic principles, namely the general principles 
of quantum theory and the lessons from general relativity: that physics is diffeomorphism 
invariant. This means that the field describing the gravitational interaction, and the geom- 
etry of spacetime is fully dynamical and interacting with the rest of the fields present; when 
one considers the description of the full gravity-matter system, this better be a background 
independent one. The fact that LQG is based in general principles of quantum mechanics 
means only that one is looking for a description based on the language of Quantum Me- 
chanics: states are elements on a Hilbert space (well defined, one expects), observables will 
be Hermitian operators thereon, etc. This does not mean that one should use all that is 
already known about quantizing fields. Quite on the contrary, the tools needed to construct 
a background independent quantization (certainly not like the quantization we know), are 
rather new. 

Another important feature of LQG is that it is the most serious attempt to perform a 
full non-perturbative quantization of the gravitational field by itself. It is an attempt to 
answer the following question: can we quantize the gravitational degrees of freedom without 
considering matter on the first place? Since LQG aims at being a physical theory, which 
means it better be falsifiable, one expects to answer that question unambiguously, whenever 
one has the theory fully developed. This is one of the main present challenges of the theory, 
namely to produce predictions that can be tested experimentally. On the other hand, one 
can consider situations in which the full knowledge of the quantum theory is not needed 
in order to describe a particular physical situation. This is precisely the case of black hole 
horizons and their entropy. 

Black hole entropy is one of the most intriguing constructs of modern theoretical physics. 
On the one hand, it has a correspondence with the black hole horizon area through the 
laws of (classical) black hole mechanics [3]. On the other hand it is assumed to have a 
quantum statistical origin given that the proper identification between entropy and area 
S = A/A£p came only after an analysis of quantum fields on a fixed background [ij . One 
of the most crucial test that a candidate quantum theory of gravity must pass then is to 
provide a mechanism to account for the microscopic degrees of freedom of black holes. It 
is not unfair to say that at the moment we have only two candidates for quantum gravity 
that have offered such an explanation: string/brane theory [B| and loop quantum gravity 
The LQG formalism uses as starting point isolated horizon (IH) boundary conditions at the 
classical level, where the interior of the BH is excluded from the region under consideration. 
In this sense, the description is somewhat effective, since part of the information about the 
interior is encoded in the boundary conditions at the horizon that in the quantum theory 
get promoted to a condition that horizon states must satisfy. There is also an important 
issue regarding this formalism. Loop quantum gravity possesses a one parameter family of 
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inequivalent representations of the classical theory labelled by a real number 7, the so called 
Barbero-Immirzi (BI) parameter (it is the analogue of the 9 ambiguity in QCD 0]). It turns 
out that the BH entropy calculation provides a linear relation between entropy and area for 
very large black holes (in Planck units) as, 

S = XA( 1 ), 

where the parameter A is independent of 7 and depends only in the counting. We have put 
the 7 dependence in the Area, since the parameter appears explicitly in the area spectrum. 
The strategy that has been followed within the LQG community is to regard the Bekenstein- 
Hawking entropy S = A/ 4 as a requirement that large black holes should satisfy. This fixes 
uniquely the value of 7 = 70 once and for all, by looking at the asymptotic behavior, provided 
that one has the 'correct counting' that provides the right value for 7 . The analytic counting 
hasprovided also an expression for the first correction term that turns out to be logarithmic 

The purpose of the present article is to provide an introduction to the main ideas behind 
the loop quantum gravity program and to one of its main applications, namely the compu- 
tation of black hole entropy. In this regard, we shall describe some recent results that have 
been performed for Planck size black holes and that complement in a precise way the ana- 
lytical computations. In particular, as we shall show, even when the black holes considered 
are outside the original domain of applicability of the formalism, one can still learn from 
these considerations. 

The structure of the paper is as follows. In Sec. [IT] we present some preliminaries, such as 
the standard geometrodynamical variables for canonical gravity, the passage to connections 
and triads and the choice of classical observables to be quantized. In Sec. ITTTJwe describe the 
loop quantum geometry formalism, including some relevant geometric operators. Sec. IIVI is 
devoted to the formalism of quantum isolated horizons. We recall the classical formalism 
and the basic steps to the quantization of the horizon theory. State counting and black hole 
entropy is the subject of Sec. [V] We and with a discussion in Sec. [VTJ 



II. PRELIMINARIES 

A. Geometrodynamics 

The first step is to introduce the basic classical variables of the theory. Since the theory is 
described by a Hamiltonian formalism, this means that the 4-dim spacetime M of the form 
M = £ x R, where S is a 3-dimensional manifold. The first thing to do is to start with 
the geometrodynamical phase space T g of Riemannian metrics q a b on £ and their canonical 
momenta 7t ab (related to the extrinsic curvature K a \, of S into M by n ab = ^fq [K ab — |g afe K), 
with q = det(q a b) and K = q ab K a b). Recall that they satisfy, 

{n ab (x),q cd (y)} = 2k 5 a [c 5% 6 3 (x, y) ; {q ab (x) , q cd (y)} = R^), ^(y)} = (1) 

General Relativity in these geometrodynamical variables is a theory with constraints, which 
means that the canonical variables {q a bi ^ ab ) do not take arbitrary values but must satisfy 
four constraints: 

H b = D a (fc ab ) « and, H = y/q [i? (3) + ^(iff 2 - n ab n ab )} « (2) 
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The first set of constraints are known as the vector constraint and what they generate (its 
gauge orbit) are spatial diffeomorphisms on E. The other constraint, the scalar constraint 
(or super-Hamiltonian) generates "time reparametrizations" . We start with 12 degrees of 
freedom, minus 4 constraints means that the constraint surface has 8 dimensions (per point) 
minus the four gauge orbits generated by the constraints giving the four phase space degrees 
of freedom, which corresponds to the two polarizations of the gravitational field. 



B. Connection Dynamics 

In order to arrive at the connection formulation, we need first to enlarge the phase space 
T g by considering not metrics qw, but the co-triads e\ that define the metric by, 

<lab = e' a e{ 5ij (3) 

where i,j = 1,2,3 are internal labels for the frames. These represent 9 variables instead of 
the 6 defining the metric q a b, so we have introduced more variables, but at the same time a 
new symmetry in the theory, namely the SO (3) rotations in the triads. Recall that a triad e l a 
and a rotated triad e'^(x) = U l j(x) e J a (x) define the same metric q a b{x), with U l j(x) G SO (3) 
a local rotation. In order to account for the extra symmetry, there will be extra constraints 
(first class) that will get rid of the extra degrees of freedom introduced. Let us now introduce 
the densitized triad as follows: 

Et = \e^ c e{e\ (4) 



where rf is the naturally defined levi-civita density one antisymmetric object. Note that 
EfEffi* = qq ab . 

Let us now consider the canonical variables. It turns out that the canonical momenta to 
the densitized triad is closely related to the extrinsic curvature of the metric, 

K= J— ^E b 3 K ab (5) 
'det(E) 



For details see [14J. Once one has enlarged the phase space from the pairs (q a b,^ ab ) to 
(£"", K J b ), the next step is to perform the canonical transformation to go to the Ashtekar — 
Barbero variables. First we need to introduce the so called spin connection T l a , the one 
defined by the derivative operator that annihilates the triad e l a (in complete analogy to the 
Christoffel symbol that defined the covariant derivative D a killing the metric). It can be 
inverted from the form, 

d [a el ] + e^ k Tie k b =0 (6) 

This can be seen as an extension of the covariant derivative to objects with mixed indices. 
The key to the definition of the new variables is to combine these two objects, namely the 
spin connection T with the object K l a (a tensorial object), to produce a new connection 

~<A\:=r a + 1 Kl (7) 

This is the Ashtekar- Barbero Connection. Similarly, the other conjugate variable will be the 
rescaled triad, 

= E«h (8) 
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Now, the pair ( 7 A^, 7 .E?) will coordinatize the new phase space T 7 . We have emphasized 
the parameter 7 since this labels a one parameter family of different classically equivalent 
theories, one for each value of 7. The real and positive parameter 7 is known as the Barbero- 
Immirzi parameter In terms of these new variables, the canonical Poisson brackets 

are given by, 

VK{x)^E){y)} = K5 b a 5)5\x,y). (9) 

and, 

{rAi(x)M(v)} = VEt{x)^E){y)} = (10) 

The new constraint that arises because of the introduction of new degrees of freedom takes 
a very simple form, 

G l = V a Et^0 (11) 

that is, it has the structure of Gauss' law in Yang-Mills theory. We have denoted by D the 
covariant defined by the connection y A l a , such that V a E1 = d a E^ + e^* J A 3 a E^. The vector 
and scalar constraints now take the form, 

K = ^A fe -(l + 7 2 )^G^0 (12) 

where F l ab = <9 a 7 A£ — <9f, 7 A' a + e l jk 7 A£ 7 Aj; is the curvature of the connection 7 A£. The other 
constraint is, 

EfE. 



S 
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det(E) 



e^ k F k ab -2(l + 1 2 )KlK^ ^0 (13) 



The next step is to consider the right choice of variables, now seen as functions of the phase 
space T 7 that are preferred for the non-perturbative quantization we are seeking. As we 
shall see, the guiding principle will be that the functions (defined by an appropriate choice 
of smearing functions) will be those that can be defined without the need of a background 
structure, i.e. a metric on E. 



C. Holonomies and Fluxes 

Since the theory possesses these constraints, the strategy to be followed is to quantize first 
and then to impose the set of constraints as operators on a Hilbert space. This is known 
as the Kinematical Hilbert Space 7^km- O ne of the main achievements of LQG is that this 
space has been rigourously defined. 

Let us start by considering the connection A a . The most natural object one can construct 
from a connection is a holonomy h a (A) along a loop a. This is an element of the gauge 
group G = SU{2) and is denoted by, 

h a (A) =Pexp (j(A a d/) (14) 

The path-order exponential of the connection. Note that for notational simplicity we 
have omitted the 'lie-algebra indices'. From the holonomy, it is immediate to con- 
struct a gauge invariant function by taking the trace arriving then at the Wilson loop 
T[a] := ±Tr P exp (£ A a ds a ). 



5 



In recent years the emphasis has shifted from loops to consider instead closed graphs T, 
that consist of N edges ej (I = 1,2, . . . , N), and M vertices i> M , with the restriction that 
there are no edges with 'loose ends'. Given a graph T, one can consider the parallel transport 
along the edges e/, the end result is an element of the gauge group gj = h(ej) e G for each 
such edge. One can then think of the connection A l a as a map from graphs to iV-copies of 
the gauge group: A\ : T — > G N . Furthermore, one can think of Ar as the configuration 
space for the graph T, that is homeomorphic to G N . 

What we are doing at the moment is to construct relevant configuration functions. In 
particular, what we need is to consider generalizations of the Wilson loops T[ct] defined 
previously, making use of the graphs and the space Ar- Every graph T can be decomposed 
into independent loops on and the corresponding Wilson loops T[an] are a particular example 
of functions defined over Ar- What we shall consider as a generalization of the Wilson loop 
are all possible functions defined over Ar- Thus, a function c : G N — > C defines a cylindrical 
function Cr of the connection A as, 

Cr :=c(h(e 1 ),h(e 2 ),...,h(e N )) (15) 

By considering all possible functions c and all possible embedded graphs T, we generate the 
algebra of functions known as Cyl (it is closely related to the holonomy algebra, and it can 
be converted into a C*-algebra Cyl, by suitable completion). 

Let us now discuss why this choice of configuration functions is compatible with the basic 
guiding principles for the quantization we are building up, namely diffeomorphism invariance 
and background independence. Background independence is clear since there is no need for 
a background metric to define the holonomies. Diffeomorphism invariance is a bit more 
subtle. Clearly, when one applies a diffeomorphism <fi : E — > S, the holonomies transform in 
a covariant way 

0* • h(ei) = h^' 1 ■ e/ ) , (16) 

that is, the diffeomorphism acts by moving the edge (or loop). How can we then end up with 
a diffeo-invariant quantum theory? The strategy in LQG is to look for a diffeomorphism 
invariant representation of the diffeo-covariant configuration functions. As we shall see later, 
this has indeed been possible and in a sense represents the present 'success' of the approach. 

Let us now consider the functions depending of the momenta that will be fundamental 
in the (loop) quantization. The basic idea is again to look for functions that are defined 
in a background independent way, that are natural from the view point of the geometric 
character of the object (1-form, 2-form, etc), and that transform covariantly with respect to 
the gauge invariances of the theory. Just as the the connection A\ can be identified with a 
one form that could be integrated along a one-dimensional object, one would like analyze the 
geometric character of the densitised triad E in order to naturally define a smeared object. 
Recall that the momentum is a density-one vector field on E, Ef with values in the dual of 
the lie-algebra su(2). In terms of its tensorial character, it is naturally dual to a (lie-algebra 
valued) two form, 

Eabi '■— \ r)abcE1 (17) 

where rj a b c is the naturally defined Levi-Civita symbol. It is now obvious that the momenta 
is crying to be integrated over a two-surface S. It is now easy to define the objects 

E[S,f] := [ E abt fdS ab , (18) 
Js 
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where f l is a lie-algebra valued smearing function on S. This 'Electric flux' variable does 
not need a background metric to be defined, and it transforms again covariantly as was the 
case of the holonomies. The algebra generated by holonomies and flux variables is known 
as the Holonomy-Flux algebra TiT. 

Perhaps the main reason why this Holonomy-Flux algebra TiT is interesting, is the way 
in which the basic generators interact, when considering the classical (Poisson) lie-bracket. 
First, given that the configuration functions depend only on the connection and the con- 
nections Poisson-commute, one expects that {T[a], T[/?]} = for any loops a and /3. The 
most interesting poisson bracket one is interested in is the one between a configuration and 
a momenta variable, 

{T[a],E[S, /]} = k £ />) S \v) Tr fa H<*)) (19) 
m 

where the sum is over the vertices v and S\ v ) = ±1 is something like the intersection 
number between the loops a and the surface S at point v. The sum is over all intersection 
of the loop a and the surface S. The most important property of the Poisson Bracket is 
that it is completely topological. This has to be so if we want to have a fully background 
independent classical algebra for the quantization. 

A remark is in order. The value of the constant t\ v depends not only on the relative 
orientation of the tangent vector of the loop a with respect to the orientation of S and S, 
but also on a further decomposition of the loop into edges, and whether they are 'incoming' 
or 'outgoing' to the vertex v. The end result is that is we have, for simple intersections, that 
the number t\ v becomes insensitive to the 'orientation'. This is different to the U(l) case 
where the final result is the intersection number. For details see [l2| • 

Let us now consider the slightly more involved case of a cylindrical function Cy that is 
defined over a graph T with edges ej, intersecting the surface S at points p. We have then, 

{C„ E[S, /]} = £££ <h) f*(p) X) p ■ c (20) 
p i P 

where the sum is over the vertices p of the graph that lie on the surface S, I p are the edges 
starting or finishing in p and where X\ p ■ c is the result of the action of the i-th left (resp. 
right) invariant vector field on the I p -th copy of the group if the I p -th edge is pointing away 
from (resp. towards) the surface S. Note the structure of the right hand side. The result is 
non-zero only if the graph T used in the definition of the configuration variable Cy intersects 
the surface 5" used to smear the triad. If the two intersect, the contributions arise from the 
action of right/left invariant vector fields on the arguments of c associated with the edges 
at the intersection. 

Finally, the next bracket we should consider is between two momentum functions, namely 
{E[S, /], E[S', g]}. Just as in the case of holonomies, these functions depend only on one 
of the canonical variables, namely the triad E. One should then expect that their Poisson 
bracket vanishes. Surprisingly, this is not the case and one has to appropriately define the 
correct algebraic structure 1 . 



The end result is that one should not regard E[S, f] as phase space functions subject to the ordinary 
Poisson bracket relations, but rather should be viewed as arising from vector fields X a on A. The non- 
trivial bracket is then due to the non-commutative nature of the corresponding vector fields. This was 
shown in 1121] where details can be found 
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We have arrived then to the basic variables that will be used in the quantization in order 
to arrive at LQG. They are given by, 

h(ei) Configuration function (21) 

and 

E[S, f] Momentum function, (22) 

subject to the basic Poisson bracket relations given by Eqs. ffTU]) and (12"U1) . In the next section 
we shall take the Holonomy-Flux algebra 7iJF as the starting point for the quantization. 



III. LOOP QUANTUM GEOMETRY 

Let us now consider the particular representation that defines LQG. As we have discussed 
before, the basic observables are represented as operators acting on wave functions \&-r(-<4) £ 
Tikin as follows: 

h{ei) ■ = (h( ei ) vl/) (A) (23) 

and 

E[S, /] • ¥ T (2) = ih {*r, 2 E[S, /]} = * 8tt | ]T ]T K ^ K ' * ( 24 ) 

v i P 

where i\ = G h, the Planck area is giving us the scale of the theory (recall that the Immirzi 
parameter 7 does not appear in the basic Poisson bracket, and should therefore not play any 
role in the quantum representation). Here we have assumed that a cylindrical function \Er T 
on a graph T is an element of the Kinematical Hilbert space (which we haven't defined yet!). 
This implies one of the most important assumption in the loop quantization prescription, 
namely, that objects such as holonomies and Wilson loops that are smeared in one dimension 
are well defined operators on the quantum theory 2 . 

The basic idea for the construction of both the Hilbert space 7ikm (with its measure) and 
the quantum configuration space A, is to consider all possible graphs on S. For any given 
graph T, we have a configuration space Ar = (SU(2)) N , which is n-copies of the (compact) 
gauge group SU(2). Now, it turns out that there is a preferred (normalized) measure on 
any compact semi- simple Lie group that is left and right invariant. It is known as the Haar 
measure /ih on the group. We can thus endow Ar with a measure fly that is defined by using 
the Haar measure on all copies of the group. Given this measure on Ar, we can consider 
square integrable functions thereon and with them the graph-T Hilbert space Tir, which is 
of the form: 

H r = L 2 (A r ,dfi T ) (25) 

If we were working with a unique, fixed graph T (which we are not), we would be in the 
case of a lattice gauge theory on an irregular lattice. The main difference between that 
situation and LQG is that, in the latter case, one is considering all graphs on E, and one 
has a family of configurations spaces {Ay /T a graph in S}, and a family of Hilbert spaces 
{Hr /T a graph inS}. 



2 this has to be contrasted with the ordinary Fock representation where such objects do not give raise to 
well defined operators on Fock space. This implies that the loop quantum theory is qualitatively different 
from the standard quantization of gauge fields. 
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The quantum configuration space A is the configuration space for the "largest graph"; 
and similarly, the kinematical Hilbert space 7ikin is the largest space containing all Hilbert 
spaces in {Hy /T a graph in S}. The Ashtekar-Lewandowski measure /xal on ^kin is then 
the measure whose projection to any Ay yields the corresponding Haar measure /iy. The 
resulting Hilbert space can thus be written as 

ftkin = L 2 (A, d/iAL) 

The cylindrical functions ^y £ Cyl belong to the Hilbert space of the theory. 

Let us then recall what is the structure of simple states in the theory. The vacuum or 
'ground state' |0) is given by the unit function. One can then create excitations by acting 
via multiplication with holonomies or Wilson loops. The resulting state |a) = T[a] ■ |0) is an 
excitation of the geometry but only along the one dimensional loop a. Since the excitations 
are one dimensional, the geometry is sometimes said to be polymer like. In order to obtain a 
geometry that resembles a three dimensional continuum one needs a huge number of edges 
(10 68 ) and vertices. 

A. A choice of basis: Spin Networks 

The purpose of this part is to provide a useful decomposition of the Hilbert space Tly, for 
all graphs. From our previous discussion we know that the Hilbert space Hy is the Cauchy 
completion of the space of cylinder functions on T, Cyl T with respect to the norm induced 
by the Haar measure on the graph configuration space Ay = (SU(2)) N . Thus, what we are 
looking for is a convenient basis for functions Fy(A) of the form 

Fy(A) :=f(h( ei ),h( e2 ),...,h(e N )) (26) 

Let us for a moment consider just one edge, say e^. What we need to do is to be able 
to decompose any function F on G (in this case we only have one copy of the group), in a 
suitable basis. 

In the case of the group G = SU(2), there is a decomposition of a function f(g) of the 
group (g e G). It reads, 

f(g) = E v'7i7~TT./7'"' Hw (g) (27) 

3 

where, 

JT' = VTO + l) / d(i H IW (g- 1 ) f(g) (28) 
Jg 

j 

and is the equivalent of the Fourier component. The functions tl mm i (g) play the role of 
the Fourier basis. In this case these are unitary representation of the group, and the label 
j labels the irreducible representations. In the SU(2) case with the interpretation of spin, 
these represent the spin-j representations of the group. In our case, we will continue to use 
that terminology (spin) even when the interpretation is somewhat different. 

Given a cylindrical function \I>x[^4] = i/>(h(ei), hfa), ■ ■ ■ , h(ejsf)), we can then write an 
expansion for it as, 

^ r [A] = iP(h(e 1 ),h(e 2 ),...,h(e N )) 

= E ff^ ,B1 ''^ AnxWei))---^^)), (29) 

31— 3 N 
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where 44nn(9) = ylU + 1) n m n (g) is the normalized function satisfying 



dfl H (j) J m n(g) (f> 3 m > n '{g) = 8j,j>8 m ,m'$n,r, 



G 

The expansion coefficients can be obtained by projecting the state \^y), 

J 31—3N \ rmim "m N n N 1^1/ V 1 -'"/ 

This implies that the products of components of irreducible representations n^Li 0m;nJM e «)] 
associated with the N edges ej G T, for all values of spins j and for — j < m,n < j and for 
any graph T, is a complete orthonormal basis for Hki n - We can the write, 

H T = ®j Hy d (31) 

where the Hilbert space Ti, a j for a single loop a, and a label j is the familiar (2? + 1) 
dimensional Hilbert space of a particle of 'spin j\ For a complete treatment see [14| . 

In the case of geometry with group SU(2), the graphs with labelling ji = j are known as 
spin networks. As the reader might have noticed, in the geometry case there are more labels 
than the spins for the edges. Normally these are associated to vertices and are known as 
intertwiners. This means that the Hilbert spaces T^tj is finite dimensional. Its dimension 
being a measure of the extra freedom contained in the intertwiners. One could then introduce 
further labelling 1 for the graph, so we can decompose the Hilbert space as 

H r = ® j H T j = Wtj,i (32) 

where now the spaces TtrAi are one-dimensional. For more details see [3, Q and 0- With 
this convenient basis it is simple to consider geometrical operator. The most important one 
in the study of black holes is given by the flux and the area operators that we consider next. 



B. Flux and Area operators 

The operators E[S, f] corresponding to the electric flux observables, are in a sense the 
basic building blocks for constructing the quantum geometry. We have seen in SecJIIII the 
action of this operators on cylindrical functions, 

E[S, /] • vD T (Z) = -i h {M> T , 2 E[S, /]} = -i Stt | E K ^ f<P) X l • ^ ( 33 ) 

p i P 

Here the first sum is over the intersections of the surface S with the graph T, and the 
second sum is over all possible edges I p that have the vertex v p (in the intersection of S and 
the graph) as initial of final point. In the simplest case of a loop a, there are only simple 
intersections (meaning that there are two edges for each vertex), and in the simplest case 
of only one intersection between S and a we have one term in the first sum and two terms 
in the second (due to the fact that the loop a is seen as having a vertex at the intersection 
point). In this simplified case we have 

E[S, /] • V a (A) = -i Svr 4 f{p) X% ■ V> (34) 
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Note that the action of the operator is to 'project' the angular momentum in the direction 
given by /* (in the internal space associated with the Lie algebra). As we shall see, this 
operator is in a sense fundamental the fundamental entity for constructing (gauge invariant) 
geometrical operators. For this, let us rewrite the action of the flux operator (I3"3"|) . dividing 
the edges that are above the surface S, as 'up' edges, and those that lie under the surface 
as 'down' edges. 



f-2 



where the sum is over the vertices at the intersection of the graph and the surface, and 
where the 'up' operator J?,^ = Jf ' 6l + Jf ,62 + ■ ■ • + Jf ,e " is the sum over all the up edges and 

the down operator Jf,^ is the corresponding one for the down edges. 

The second simplest operator that can be constructed representing geometrical quantities 
of interest is the area operator, associated to surfaces S. The reason behind this is again 
the fact that the densitized triad is dual to a two form that is naturally integrated along a 
surface. The difference between the classical expression for the area and the flux variable 
is the fact that the area is a gauge invariant quantity. Let us first recall what the classical 
expression for the area function is, and then we will outline the regularization procedure to 
arrive at a well defined operator on the Hilbert space. The area A[S] of a surface S is given 
by A[S] = J s d 2 x yh, where h is the determinant of the induced metric h ab on S. When the 
surface S can be parametrized by setting, say, a; 3 = 0, then the expression for the area in 
terms of the densitized triad takes a simple form: 



A[S] = 7 J d 2 x ^Ef E] fcv (36) 

where = 5^ is the Killing-Cartan metric on the Lie algebra, and 7 is the Barbero-Immirzi 
parameter (recall that the canonical conjugate to A is ^Ef = Ef/j). Note that the functions 
is again smeared in two dimensions and that the quantity inside the square root is very 
much a square of the (local) flux. One expects from the experience with the flux operator, 
that the resulting operator will be a sum over the intersecting points p, so one should focus 
the attention to the vertex operator 



k l] 



(TP _ TP )( TP _ TP > 

with this, the area operator takes the form, 

A[S]=8n 1 f P J2V^'p ( 3? ) 

p 

We can now combine both the form of the vertex operator with Gauss' law (Jf^ + Jf^))'^ = ® 
to arrive at, 

l(^«)--? W )r=|2(4,)| 2 (38) 

where we are assuming that there are no tangential edges. The operator J p , » is an angular 
momentum operator, and therefore its square has eigenvalues equal to j u (j u + 1) where j u is 
the label for the total 'up' angular momentum. We can then write the form of the operator 



i[£].jV(T,i) = 87r 7 4E M«)\ 2 -W>3) (39) 
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With these conventions, in the case of simple intersections between the graph T and the 
surface S, the operator takes the well known form: 

v&V 

when acting on a spin network J\l "(7 ', j) defined over T and with labels j on the edges (we 
have not used a label for the intertwiners) . As we shall se when we consider the quantum 
theory of isolated horizons, the two operators considered here will play an important role 
not only in the geometry of the horizon but in the entropy counting. 



IV. QUANTUM ISOLATED HORIZONS 

Let us focus on the sector of the theory consisting of space-times which admit a type I 
isolated horizon A with a fixed cir<3cl Qjq clS the internal boundary. Then £ is asymptotically 
flat and has an internal boundary S, topologically a 2-sphere, the intersection of £ with A. 
The type I isolated horizon boundary conditions require that i) A be null, ii) Non-expanding, 
iii) The field equations be satisfied there and iv) the intrinsic geometry on A be left invariant 



by the null vecto i a generating A. For details see [15 . 

Introduce on S an internal, unit, radial vector field r % (i.e. any isomorphism from the unit 
2-sphere in the Lie algebra of SU(2) to S). Then it turns out that the intrinsic geometry of S 
is completely determined by the pull-back A^ri =:W to S of the (internal-radial component 
of the) connection A 1 on £ jl5| . Furthermore, W is in fact a spin-connection intrinsic to the 
2-sphere S. Finally, the fact that S is (the intersection of £ with) a type I isolated horizon 
is captured in a relation between the two canonically conjugate fields: 

F : = dW = - — £>;. (40) 
a Q 

where £ 4 is the pull-back to S of the 2-forms T, l ab = r] a b c E lc on £. (Throughout, = will stand 
for equality restricted to A.) Thus, because of the isolated horizon boundary conditions, 
fields which would otherwise be independent are now related. In particular, the pull-backs to 
S of the canonically conjugate fields A 1 , £ l are completely determined by the U(l) connection 
W. 

In absence of an internal boundary, the symplectic structure is given just by a volume 



integral [15j. Jn presence of the internal boundary under consideration, it now acquires a 
surface term 



n(Si,6 2 ) = [ Tt(MA<5 2 £-MA<Ji£) + — / StWA S 2 W , (41) 

8vrG IJ M 77T Js 

where 5 = (5 A, 5T,) denotes tangent vectors to the phase space T. Since W is essentially the 
only 'free data' on the horizon, it is not surprising that the surface term of the symplectic 
structure is expressible entirely in terms of W. However, it is interesting that the new 
surface term is precisely the symplectic structure of a U(l)-Chern Simons theory. The 
symplectic structures of the Maxwell, Yang-Mills, scalar and dilatonic fields do not acquire 
surface terms and, because of minimal coupling, do not feature in the gravitational symplectic 
structure either. Conceptually, this is an important point: this, in essence, is the reason 
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why the black hole entropy depends just on the horizon area and not, in addition, on the 
matter charges p]. 

One can systematically 'quantize' this sector of the phase space 0. We can focus only 
on the gravitational field since the matter fields do not play a significant role. One begins 
with a Kinematic Hilbert space TC = TCy ® TCs where TCy is the Hilbert space of states 
in the bulk as described before and TCs is the Hilbert space of surface states. Expression 
( [H]) of the symplectic structure implies that TCs should be the Hilbert space of states of 
a Chern-Simons theory on the punctured S, where the 'level', or the coupling constant, is 
given by: 

* = (42) 



A pre-quantization consistency requirement is that k be an integer 

Our next task is to encode in the quantum theory the fact that A is a type I horizon with 
area a . This is done by imposing the horizon boundary condition as an operator equation: 

(1®#)*= - ( — (E-r)® 1 ) *, (43) 



on admissible states ^ in TC. Now, a general solution to (143]) can be expanded out in a basis: 

is an eigenvector of the 'triad operator' — (2TT'~f/a ) (S • 

r)(x) on TCy and tyg is an eigenvector of the 'curvature operator' F(x) on TCs with same 
eigenvalues. Thus, the theory is non-trivial only if a sufficiently large number of eigenvalues 
of the two operators coincide. Since the two operators act on entirely different Hilbert spaces 
and are introduced quite independently of one another, this is a very non-trivial requirement. 

Now, in the bulk Hilbert space TCy, the eigenvalues of the 'triad operator' are given by 
3: 




m I 8^{x,p I )r] ab j , (44) 

where mi are half integers and 7]^ is the natural, metric independent Levi-Civita density 
on 5* and pi are points on S at which the polymer excitations of the bulk geometry in the 
state puncture S. A completely independent calculation involving just the surface 
Hilbert space TCs, yields the following eigenvalues of F(x): 

\ njS^pj) = 2vr^i ( 45 ) 



i 



where nj are integers modulo k. Thus, with the identification —2mj = nj mod k, the two 
sets of eigenvalues match exactly. There is a further requirement or constraint that the 
numbers mi should satisfy, namely, 



^m / = (46) 



This constraint is sometimes referred to as the projection constraint, given that the 'total 
projection of the angular momentum' is zero. Note that in the Chern-Simons theory the 
eigenvalues of F(x) are dictated by the 'level' k and the isolated horizon boundary conditions 
tie it to the area parameter a Q just in the way required to obtain a coherent description of 
the geometry of the quantum horizon. 
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In the classical theory, the parameter a Q in the expression of the surface term of the 
symplectic structure (jUJ) and in the boundary condition ([4*0]) is the horizon area. However 
in the quantum theory, a Q has simply been a parameter so far; we have not tied it to the 
physical area of the horizon. Therefore, in the entropy calculation, to capture the intended 
physical situation, one constructs a suitable 'micro-canonical' ensemble. This leads to the 
last essential technical step. 

Let us begin by recalling that, in quantum geometry, the area eigenvalues are given by, 

87r 7 4iZ)vS57+l)- 
i 

We can therefore construct a micro-canonical ensemble by considering only that sub-space 
of the volume theory which, at the horizon, satisfies: 

a - 5 < Stt^I, J2 VJiWi + 1) < a o + 5 (47) 

where I ranges over the number of punctures, jj is the spin label associated with the puncture 
Pi 0- 3 

Quantum Einstein's equations can be imposed as follows. The implementation of the 
Gauss and the diffeomorphism constraints is the same as in 0. The first says that the 
'total' state in TC be invariant under the SU(2) gauge rotations of triads and, as in [6], this 
condition is automatically met when the state satisfies the quantum boundary condition 
(|43p . The second constraint says that two states are physically the same if they are related 
by a diffeomorphism. The detailed implementation of this condition is rather subtle because 
an extra structure is needed in the construction of the surface Hilbert space and the effect of 
diffeomorphisms on this structure has to be handled carefully [6J. However, the final result 
is rather simple: For surface states, what matters is only the number of punctures; their 
location is irrelevant. The last quantum constraint is the Hamiltonian one. In the classical 
theory, the constraint is differentiable on the phase space only if the lapse goes to zero on 
the boundary. Therefore, this constraint restricts only the volume states. However, there 
is an indirect restriction on surface states which arises as follows. Consider a set (pi,ji) 
with I = 1, 2, . . . N consisting of N punctures pi and half-integers jj, real, satisfying (1471) . 
We will refer to this set as 'surface data'. Suppose there exists a bulk state satisfying the 
Hamiltonian constraint which is compatible with this 'surface data'. Then, we can find 
compatible surface states such that the resulting states in the total Hilbert space Ti lie in 
our ensemble. The space S , ( Pr j r ) of these surface states is determined entirely by the surface 
data. In our state counting, we include the number Af( Pl ,jj) of these surface states, subject 
however, to the projection constraint that is purely intrinsic to the horizon. 4 



3 The appearance of the parameter 6 is standard in statistical mechanics. It has to be much smaller than 
the macroscopic parameters of the system but larger than level spacings in the spectrum of the operator 
under consideration. 

4 Note that there may be a large number -possibly infinite- of bulk states which are compatible with a 
given 'surface data' in this sense. This number does not matter because the bulk states are 'traced out' in 
calculating the entropy of the horizon. What matters for the entropy calculation is only the dimensionality 
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V. BLACK HOLE ENTROPY 



In this section we shall deal with the issue of entropy counting. We have started with a type I 
isolated horizon of area a Q (in vacuum this is the only multipole defining the horizon), and we 
have quantized the theory and arrived to a Hilbert space as described before. The question 
now is: How many microstates correspond to the given macrostate, defined uniquely by a D ? 
Let us now pose the condition that the states in SfajA should satisfy: 

• They belong to the physical Hilbert space on the surface TCs- 

• The condition fH7|) is satisfied. 

• The quantum boundary condition (PI51) is satisfied. 

• The projection constraint fH6|) is satisfied. 

In terms of a concrete counting the problem is posed a follows: We shall consider the 
lists (p/,j/,m/) corresponding to the allowed punctures, spins of the piercing edges, and 
'projected angular momentum' labels, respectively. 

The task is then to count these states and find A/ r ( P/J/im/ ). The entropy will be then, 

S BH := HM(piji,mi)) (48) 



This problem was systematically addressed in [6| in the approximation of large horizon area 
a a . Unfortunately, the number of such states was underestimated in |6[. 5 In [9j the counting 
was completed and it was shown in detail that, for large black holes (in Planck units), the 
entropy behaves as: 

Sbu = j ~ - In A, 

provided the Barbero-Immirzi parameter 7 is chosen to coincide with the value 70, that has 
to satisfy [§]: 

1 = ^(2 h + 1) exp (-27r 7 ovSi7Tl)) • (49) 
1 

The solution to this transcendental equation is approximately 70 = 0.27398. . . (see j§, [TtJ 
for details). 

Here we shall perform the counting in a different regime, namely for small black holes 
in the Planck regime [3]. Thus we shall perform no approximations as in the previous 
results. Thus our results are complementary to the analytic treatments. On the other hand 
our counting will be exact, since the computer algorithm is designed to count all states 
allowed. Counting configurations for large values of the area (or mass) is extremely difficult 
for the simple reason that the number of states scales exponentially. Thus, for the computing 
power at our disposal, we have been able to compute states up to a value of area of about 
a Q = 550 Ip (recall that the minimum area gap for a spin 1/2 edge is about a m \ n m Qlp, 
so the number of punctures on the horizon is below 100). At this point the number of 



5 The under-counting was noticed in and a new counting was there proposed and carried out in 
However the choice of relevant states to be counted there is slightly different from our case. Details of the 
comparison between two methods will be reported elsewhere [17| . 
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FIG. 1: The entropy as a function of area is shown, where the projection constraint has not been 
imposed. The BI is taken as 7 = 0.274. 

states exceeds 2.8 x 10 58 . In terms of Planck masses, the largest value we have calculated is 
M = 3.3 Mp. When the projection constraint is introduced, the upper mass we can calculate 
is much smaller, given the computational complexity of implementing the condition. In this 
case, the maximum mass is about 1.7 Mp. 

It is important to describe briefly what the program for counting does. What we are 
using is what it is known, within combinatorial problems, as a brute force algorithm. This 
is, we are simply asking the computer to perform all possible combinations of the labels 
we need to consider, attending to the distinguishability -indistinguishability criteria that 



are relevant [6|, |17|, and to select (count) only those that satisfy the conditions needed to 
be considered as permissible combinations, i.e., the area condition and the spin projection 
constraint. An algorithm of this kind has an important disadvantage: it is obviously not the 
most optimized way of counting and the running time increases rapidly as we go to little 
higher areas. This is currently the main limitation of our algorithm. But, on the other 
hand, this algorithm presents a very important advantage, and this is the reason why we 
are using it: its explicit counting guarantees us that, if the labels considered are correct, the 
result must be the right one, as no additional assumption or approximation is being made. 
It is also important to have a clear understanding that the algorithm does not rely on any 
particular analytical counting available. That is, the program counts states as specified in 
the original formalism [6]. The computer program has three inputs: i) the classical mass M 
(or area a Q = lQir M 2 ), ii) The value of 7 and iii) The size of the interval S. 

Once these values are given, the algorithm computes the level of the horizon Chern Simons 
theory k = [a /47T7] and the maximum number of punctures possible n max = [a /47T7\/3] 
(where [■], indicates the principal integer value). At first sight we see that the two conditions 
we have to impose to permissible combinations act on different labels. The area condition 
acts over j's and the spin projection constraint over m's. This allows us to first perform 
combinations of j's and select those satisfying the area condition. After that, we can perform 
combinations of m's only for those combinations of j's with the correct area, avoiding some 
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FIG. 2: Entropy vs Area with and without the projection constraint, with 5 = 0.5. 

unnecessary work. We could also be allowed to perform the counting without imposing the 
spin projection constraint, by simply counting combinations of j's and including a multi- 
plicity factor of n/(2j/ + 1) f° r each one, accounting for all the possible combinations of 
m's compatible with each combination of j's. This would reduce considerably the running 
time of the program, as no counting over m's has to be done, and will allow us to separate 
the effects of the spin projection constraint (that, as we will see, is the responsible of a 
logarithmic correction). It is very important to notice at this point that this separation of 
the counting is completely compatible with the distinguishability criteria. 

The next step of the algorithm is to consider, in increasing order, all the possible number 
of punctures (from 1 to n max ) and in each case it considers all possible values of jj. Given 
a configuration j 2 , . . . , j n ) {n < n max ), we ask whether the quantum area eigenvalue 
A = 87T7y 'ji (jj + 1) lies within [a Q — 8, a a + 8]. If it is not, then we go to the next 
configuration. If the answer is positive, then the labels m's are considered as described 
before. That is, for each of them it is checked whether mj = is satisfied. 

In Figure [H we have plotted the entropy, as S = ln(# states) vs the area a Q , where we 
have counted all possible states without imposing the Yl m i constraint, and have chosen a 
5 = 0.5. As it can be seen, the relation is amazingly linear, even for such small values of the 
area. When we fix the BI parameter to be 7 = 70 = 0.274, and approximate the curve by a 
linear function, we find that the best fit is for a slope equal to 0.2502. 

When we include the projection constraint, the computation becomes more involved and 
we are forced to consider a smaller range of values for the area of the black holes. In Figure [21 
we plot both the entropy without the projection and with the projection, keeping the same 
5. The first thing to note is that for the computation with the constraint implemented, there 
are some large oscillations in the number of states. Fitting a straight line gives a slope of 
0.237. In order to reduce the oscillations, we increased the size of 8 to 5 = 2. The result 
is plotted in Figure [3j As can be seen the oscillations are much smaller, and the result of 
implementing the constraint is to shift the curve down (the slope is now 0.241). In order to 
compare it with the expected behavior of S = A/4 — (1/2) In A, we subtracted both curves 
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FIG. 3: Entropy vs Area with and without the projection constraint, with 5 = 2. 

of Figure El in the range a Q = [50,160], and compared the difference with a logarithmic 
function. The coefficient that gave the best fit is equal to —0.4831 (See Figure Hj). What 
can we conclude from this? While it is true that the rapidly oscillating function is far from 
the analytic curve, it is quite interesting that the oscillatory function follows a logarithmic 
curve with the "right" coefficient. It is still a challenge to understand the meaning of the 
oscillatory phase. Even when not conclusive by any means, we can say that the counting of 
states is consistent with a (n asymptotic) logarithmic correction with a coefficient equal to 
(-1/2). 

VI. DISCUSSION 

In this contribution we have considered the approach to the quantum theory of gravity 
known as loop quantum gravity. We have presented a brief introduction to the main ideas 
behind this approach and have considered one of its main applications, black hole entropy. 
We have discussed the main features in the approach to black hole entropy, in particular in 
the implementation of the isolated horizon boundary conditions to the quantum theory and 
how this conditions tell us what states can be regarded as 'black hole states' that contribute 
to the entropy of the horizon. As we have seen, the fact that there is an intrinsic discreteness 
in the quantum horizon theory and that we are ignoring (tracing out) the states in the bulk, 
is the reason why the entropy becomes finite. It is sometimes believed that the fact that 
we do get an entropy proportional to area is natural are not surprising, given that on the 
horizon, the theory under consideration is a Chern-Simons theory with punctures, and the 
entropy of a two dimensional theory should be proportional to the total volume (area in this 
case). It is important to stress that the result is not as trivial as it sounds. To begin with, 
we do not have a given Chern-Simons theory on the horizon, for any macro-state of a given 
area a Q , there are many possible microstates that can be associated with it. They do not all 
live on the same 'Chern-Simons Hilbert space'. The surface Hilbert space Tis is made out 



18 



-3,2 

50 60 70 80 90 100 110 120 130 140 150 

Area 




FIG. 4: The curves of Fig. [3] are subtracted and the difference, an oscillatory function, shown in 
the upper figure. The curve is approximated by a logarithm curve in the lower figure. 



of the tensor product of all possible Chern-Simons states compatible with the constraints 
detailed in Sec. refsec:5, which belong to different Chern-Simons states (characterized by, 
say, the total number of punctures). That the total number of states compatible with the 
(externally imposed) constraints is proportional to area is thus a rather non-trivial result. 

One might also wonder about the nature of the entropy one is associating to the horizon, 
there has been some controversy about the origin and location of the degrees of freedom 
responsible for black hole entropy (see for instance 3] for a recent discussion). It has been 
argued that the degrees of freedom lie behind the horizon, on the horizon and even on an 
asymptotic region at infinity. What is then our viewpoint on this issue? The viewpoint is 
that the IH boundary conditions implement in a consistent manner an effective description, 
as horizon data, of the degrees of freedom that might have formed the horizon. These 
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degrees of freedom, even if they are there in the physical space-time, they are inaccessible to 
an external observer. The only thing that this observer can 'see' are the degrees of freedom 
at the horizon, and these degrees of freedom are thus responsible for the entropy associated 
to the horizon. 

In the last part of this article, we have focussed our attention on some recent results 
pertaining to the counting of states for Planck size horizons. As we have shown, even when 
these black holes lie outside the original domain of validity of the isolated horizon formalism 
(tailored for large black holes), the counting of such states has shed some light on such 
important issues as the BI parameter, responsible for the asymptotic behavior, and the first 
order, logarithmic, correction to entropy. We have also found, furthermore, that there is a 
rich structure underlying the area spectrum and the number of black hole states that could 
not have been anticipated by only looking at the large area limit. In particular, it has been 
found that the apparent periodicity in the entropy vs area relation yields an approximate 
'^quantization of the entropy that makes contact with Bekenstein's heuristic considerations 
20| , and is independent on the choice of relevant states, and its associated counting. Details 



will be published elsewhere [21 
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